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(Dated: February 3, 2008) 

We present measurements of the total production rates and momentum distributions of the 
charmed baryon yl+ in e^e~ — > hadrons at a center-of-mass energy of 10.54 GeV and in T(45') 
decays. In hadronic events at 10.54 GeV, charmed hadrons are ahnost exchisively leading particles 
in e^e^ ^ cc events, allowing direct studies of c-quark fragmentation. We measure a momentum 
distribution for baryons that differs significantly from those measured previously for charmed 
mesons. Comparing with a number of models, we find none that can describe the distribution 
completely. We measure an average scaled momentum of (xp) = 0.574±0.009 and a total rate of 
A'^^^= 0.057±0.002(exp.)±0.015(BF) yljt per hadronic event, where the experimental error is much 
smaller than that due to the branching fraction into the reconstructed decay mode, pK^ix^ . In 
r(4S') decays we measure a total rate of N^^ = 0.091±0.006(exp.)±0.024(BF) per r(4S') decay, 
and find a much softer momentum distribution than expected from B decays into a plus an 
antinucleon and one to three pions. 

PACS numbers: 13.66.Bc, 13.20.He, 13.60.Rj 



I. INTRODUCTION 

The production properties of charmed baryons in e+e~ 
annihilations into cc and in decays of bottom (6) hadrons 
probe different aspects of strong interaction physics. Ex- 
periments running at and below the T{AS) resonance are 
uniquely positioned to explore each of these processes in 
detail. The CLEO experiment has made precise stud- 
ies of charmed mesons in this way [1], and the larger 
data samples available at the B factories have allowed 
improved studies of charmed mesons [2] and the first pre- 
cise studies of charmed baryons [2, 3]. 

Heavy hadrons [H) produced in e+e~ annihilations 

provide a laboratory for the study of heavy-quark {Q = 
c, b) jet fragmentation, in terms of both the relative pro- 
duction rates of hadrons with different quantum num- 
bers and their associated spectra. The latter can be 
characterized in terms of a scaled energy or momentum, 
such as Xp = P*[]/p*nax^ wlierc p*^ is the hadron mo- 
mentum in the e+e^ center-of-mass (cm.) frame and 
Pmax — \/s/4 — rn^j is the maximum momentum avail- 
able to a particle of mass rriH at a cm. energy of y^. 
For -^i ^ 2mH it has been observed [4] that the Xp distri- 
butions P{xp) for heavy hadrons peak at relatively high 
values, and that very few h hadrons are produced apart 
from those containing initial h quarks, so that one can 
probe leading 6-hadron production directly. This is also 
the case for charmed (c) hadrons when ^/s < 2m b-, where 
rriB is the mass of the lightest b meson, but above BB 
threshold a large fraction of the c hadrons are 6-hadron 
decay products. 

Since the hadronization process is intrinsically non- 
perturbative, P{xp) cannot be calculated using pertur- 
bative Quantum Chromodynamics (QCD). However, a 
high quark mass provides a convenient cut-off point and 
the distribution of the scaled momentum of the heavy 
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quark before hadronization, xq = 2pg/y^, can be cal- 
culated [5-8]. The observable P{xp) is thought to be re- 
lated by a simple convolution or hadronization model. 
Several phenomenological models of heavy-quark frag- 
mentation have been proposed [9 12]. Predictions de- 
pend on the mass of the heavy quark, with P{xp) be- 
ing much harder for b hadrons than c hadrons, and in 
some cases on the mass and quantum numbers of H. 
Hadrons containing the same heavy quark type are gen- 
erally predicted to have quite similar P(xp), although 
differences between mesons and baryons have been sug- 
gested [13, 14]. Measurements of P{xp) serve to constrain 
perturbative QCD and these model predictions. Further- 
more, measurements for a given c or 6 hadron at different 
y/s can test QCD evolution, and comparisons of c- and b- 
hadron distributions can test heavy-quark symmetry [15]. 

The inclusive &-hadron scaled energy distribution and 
its average value of 0.71 have been measured precisely [16] 
by experiments at the Z", using partial reconstruction 
techniques. However, these techniques do not distinguish 
the different types of b hadrons. The relative production 
of B~ , B^, B^, excited b mesons and b baryons have been 
measured [4, 17], but with limited precision and no sen- 
sitivity to differences in their Xp distributions. Several c 
mesons have been studied at the Z° [18], but it is diffi- 
cult to disentangle the leading charm and 6-decay con- 
tributions and neither component is measured precisely. 
Recent measurements below BB threshold [1, 2] have 
good precision over the full Xp range and show substan- 
tial differences between the pseudoscalar D and vector 
D* meson states. P{xp) has been measured below BB 
threshold for two charmed baryons, A'^ by CLEO [19] and 
Belle [2], and S"" by BaBAR [3], but with limited statistics, 
especially at low Xp. In this article we use the excellent 
particle identification of the BaBAR experiment to iso- 
late yl^ baryons (the inclusion of charge conjugate states 
is implied throughout) in a large data sample collected 
below BB threshold, at \/s = 10.54 GeV. We measure 
P{xp) precisely, and compare our results with available 
predictions and previous measurements of heavy hadrons. 

The large &-hadron masses allow many hadronic decay 
modes, a small fraction of which have been studied in de- 
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tail. The T(45) resonance provides a unique laboratory, 
in which no b baryons are produced and decays of mesons 
into baryons can be studied directly. Many c baryons 
have been observed in inclusive T(4S') decays [4] and the 
low rate of associated leptons and high rate of "wrong- 
sign" yl+ [20] suggest interesting dynamics. However 
only a few exclusive decays with c baryons have been 
observed [21-23]. Again, momentum distributions have 
been measured only for the /!+ [2, 24] and [3] with lim- 
ited precision. Here we use data collected on the T(4S') 
resonance {^/s = 10.58 GeV) and subtract the e+e~— >cc 
contribution measured at \/s = 10.54 GeV to make a pre- 
cise measurement of the momentum distribution in 
B meson decays, which we compare with a number of 
possible models. 

In section H, we describe the BaBAR detector, in par- 
ticular the particle identification capabilities essential to 
these measurements. In sections III and IV, we discuss 
the selection of Af candidates and the measurement of 
their Xp distributions, respectively. We interpret the re- 
sults for e+e~^gg events and T(4S') decays in sections V 
and VI, respectively, and summarize in section VII. 



II. APPARATUS AND EVENT SELECTION 

In this analysis, we use data samples corresponding to 
9.5 fb~^ of integrated luminosity at ^/s = 10.54 GeV 
and 81 fb'^ on the 7(45*) resonance, ^=10.58 GeV. 
The BaBAR detector is located at the PEP-II asymmetric- 
energy e+e~ collider at the Stanford Linear Accelerator 
Center and is described in detail in Ref. [25]. We use 
charged tracks measured in the five-layer, double sided 
silicon vertex tracker (SVT) and the 40-layer drift cham- 
ber (DCH). In a 1.5 T axial magnetic field, they provide 
a combined resolution on the momentum pT transverse 
to the beam axis of [ct(pt)/pt]^ = [0.0013pt]^ +0.0045^, 
where pt is measured in GeV/ c. 

Charged particle identification uses a combination 
of the energy loss (dE / dx) measured in the DCH, 
and information from the detector of internally re- 
fiected Cherenkov light (DIRC). The DCH gas is he- 
lium:isobutane 80:20 and a typical cell size is 18 mm. A 
truncated mean algorithm gives a dE/dx value for each 
track with an average resolution of 7.5%, from which we 
calculate a set of five relative likelihoods if*"^^ for the 
particle hypotheses i = e,^,7r,X, and p. Differences be- 
tween the log-likehhoods ifp^ = ln(Lpc") - \n{L^^^) 
are used as input to the particle identification algorithm. 

The DIRC comprises 144 fused silica bars that guide 
Cherenkov photons to an expansion volume filled with 
water and equipped with 10752 photomultiplier tubes. 
The fused silica refractive index is 1.473, corresponding 
to Cherenkov thresholds of 128, 458 and 867 MeV/c for 
pions, kaons and protons, respectively. A particle well 
above threshold yields 20-75 measured photons, each 
with a Cherenkov angle resolution of about 10 mrad. 
A global likelihood algorithm considers all the recon- 



structed charged tracks and detected photons in each 
event and assigns each track a set of likelihoods Lf^^'^ ■ 

The DCH provides excellent K-tt and p-K separa- 
tion for momenta in the laboratory frame below 0.5 and 
0.9 GeV/c, respectively. The DIRC provides very good 
separation for momenta above 1.0 and 1.5 GeV/c, respec- 
tively. In both detectors the separation power is lower for 
tracks with polar angles near 90° in the laboratory than 
for more forward or backward tracks. To minimize the 
systematic errors in this analysis, the identification effi- 
ciencies must not vary rapidly as a function of momen- 
tum or polar angle. We have therefore developed an al- 
gorithm that uses linear combinations of ifj"^ and Z]]^^*^ 
chosen to minimize such variations. It is described in de- 
tail in Ref. [26] and its performance for tracks used in this 
analysis is shown as a function of momentum in Fig. 1. 
The identification efficiencies are better than 99% at low 
momenta and above 90% for the majority of yl+ decay 
products. They are seen to vary smoothly with momen- 
tum, and are almost independent of polar angle except 
near 0.8 GeV/c (1.2 GeV/c) for pions and kaons (pro- 
tons), where they are as much as 10% lower for central 
tracks than for forward/backward tracks. The misiden- 
tification rates depend strongly on polar angle in the 
momentum regions 0.6-0.8, 1.1-1.3 and 2.5-3.5 GeV/c. 
They are below 5% everywhere except that the rate for 
kaons (protons) to be misidentified as pions reaches 11% 
at 0.7 (1.2) GeV/c for the most central tracks. These 
rates have negligible effects on the results. About 16% of 
the selected tracks have good DCH information but are 
outside the DIRC fiducial acceptance. These can be iden- 
tified with essentially the same efficiencies as in Fig. 1 for 
pion and kaon (proton) momenta below 0.6 (0.9) GeV/c. 

The event selection requires thrc^c or more charged 
tracks in the event, which retains any e'^e~^qq event or 
T(45) decay containing a reconstructable A^^pK^-K^ 
decay and suppresses beam- related backgrounds. We 
evaluate its performance using a number of simulations, 
each consisting of a generator for a certain type of event 
combined with a detailed simulation of the BaBM. detec- 
tor [27]. For e+e-^qq events we use the JETSET [28] 
generator and for T(4S') events we use our own gener- 
ator, EVTGEN [29], in which the T{4.S) decays into a 
BB pair, then the B and B decay using a combination 
of measured exchisive and semi-exclusive modes, and a 
&— » cW~ model tuned to the world's inclusive data. We 
study large samples of simulated two-photon, r-pair and 
radiative e- and /i-pair events, and find their contribu- 
tions to both signal and background to be negligible. 



III. A+^pK-TT+ SELECTION 

We construct yl+ candidates from charged tracks that are 
consistent with originating at the e"'"e~ interaction point 
and have good tracking and particle identification infor- 
mation. Each track must have: (i) at least 20 measured 
coordinates in the DCH; (ii) at least 5 coordinates in the 
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FIG. 1: EfRcicncics for identifying (solid symbols, left-hand vertical scale) or rnisidentifying (open symbols, riglit-liand scale) 
selected pions (circles), kaons (squares) and protons (triangles) within the DIRC acceptance as pions (left), kaons (center) or 
protons (right). They are extracted from control samples in the data as smooth functions of momentum and polar angle, and 
shown in momentum bins averaged over the polar angle range used. The error bars indicate the average uncertainty due to 
control sample statistics. 



SVT, including at least three in the direction along the 
e~ beam; (iii) a distance of closest approach to the beam 
axis below 1 mm; and (iv) a z-coordinatc at this point 
within 10 cm of the nominal interaction point. These 
criteria ensure good quality information from the DCH 
and a well measured entrance angle into the DIRC. If the 
extrapolated trajectory intersects a DIRC bar then the 
track is accepted if it is identified as a pion, kaon or pro- 
ton by the combined DCH and DIRC algorithm. If not, 
then it is accepted if it is identified as a proton (pion or 
kaon) using DCH information only and has a momentum 
below 1.2 (0.6) GeV/c. 

We consider a combination of three charged tracks as 
a candidate if the total charge is +1. one of the pos- 
itively charged tracks is identified as a proton and the 
other as a pion, and the negatively charged track is iden- 
tified as a kaon. With the appropriate particle type as- 
signed to each track, we correct their measured momenta 
for energy loss and calculate their combined four mo- 
mentum from their momenta at their points of closest 
approach to the beam axis. The distributions of invari- 
ant mass for the candidates in the on- and off-resonance 
data are shown in Fig. 2; signals of about 137,000 
and 13,000 decays, respectively, are visible over nearly 
uniform backgrounds. 

The yl+ reconstruction efficiency depends primarily on 
the momenta p and polar angles 6 of the daughter tracks 
in the laboratory frame. To reduce systematic uncer- 
tainty, we apply an efficiency correction to each candi- 
date before boosting it into the e+e" cm. frame. The 
efficiencies for reconstructing and identifying tracks from 
pions, kaons and protons are determined from large con- 



trol samples in the data as two-dimensional functions of 
{p, 9). We use these efficiencies in dedicated simulations 
of qq and T(4S') events containing a yl+ baryon that is 
decayed into pK~7T^. From these we calculate the yl+ se- 
lection efficiency e as a smooth two-dimensional function 
of (p, 9) of the yl+ . We check that the efficiency does not 
depend on other track or event variables, in particular 
that it is the same in simulated qq and T(4S') events for 
given values of {p, 9) . The resolutions on the momen- 
tum and polar angle are much smaller than the bin sizes 
used below, so we include resolution effects by defining 
the efficiency as the number of reconstructed within 
a given (p, 9) range divided by the number generated in 
that range, using ranges smaller than the relevant bin 
sizes. We test the efficiency using a number of simula- 
tions, and find biases to be below 1%. 

The efficiency varies rapidly near the edges of the de- 
tector acceptance and at very low momenta in the labora- 
tory. We make the tight fiducial requirement that 6*, the 
polar angle of the candidate in the e^e~ cm. frame, 
satisfy — 0.7< cos6'*<0.2, which reduces model depen- 
dence and rejects all candidates in regions with efficiency 
below about 5%, including those with a total laboratory 
momentum below about 0.7 GeV/c. A feature of the 
boosted cm. system is that true yl+ baryons with low 
cm. momentum p* are boosted forward and often have 
all three tracks in the detector acceptance, giving efficient 
access to the full p* range. 

We define Ak as the fraction of the A'^ in events of type 
k produced within our fiducial range — 0.7< cos^*<0.2. 
In r{AS)->BB-fA+X decays, the true cos 6** distribu- 
tion is uniform and Ay = 0.45. In cc events the angular 
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FIG. 2: Invariant mass distributions for candidates in the 
on- (black) and off-resonance (gray) data. 

distribution of the initial c-quark follows 1 -|- cos^ 0*, and 

we use the JETSET simulation to calculate the distri- 
bution for yl+ after QCD radiation and hadronization. 
Soft A'^ are produced predominantly in events with hard 
gluon radiation, which flattens the distribution consider- 
ably. The resulting value of A^c is 0.46 at p*=0, and falls 
with increasing p* toward an asymptotic value of 0.38. 

We bin candidates according to their reconstructed val- 
ues of Xp — p* /Pmaxj where p^ax is calculated for each 
event from its cm. energy and the nominal A'^ mass [4]. 
Figure 3 shows the average value in each Xp bin of the 
product Acc{p*) ■ £{p,6) for selected candidates in the 
off- resonance data. It ranges from 8% at low Xp to 19% 
at high Xp. The error bars represent the statistical un- 
certainty on the cfflciency calculation. The correspond- 
ing quantity for yl+ from T{AS) decays, {Ar ■ e{p, 6*)), is 
slightly higher at low Xp due to a small dependence of 
e on 9, and rises faster with increasing Xp since Ax is 
constant, whereas Acc decreases. We give each candidate 
a weight equal to the inverse of either Acc ■ £ or Ar ■ e, 
as specified below. The RMS deviation of the weights in 
each bin is always much smaller than the average value. 




FIG. 3: Average value of Acc ■ £, the A'^ —> pK'-K^ accep- 
tance times reconstruction efficiency, for candidates in the 
off-resonance data in bins of scaled momentum Xp. The error 
bars represent the average statistical uncertainty. 



The simulation predicts a nearly uniform background 
over the pK^ir^ mass range shown in Fig. 2. We search 
for reflections in the data by changing the particle mass 
assignments. We observe signals for D+^iir~7r+7r+ (7r+ 
misidentified as p) and Df^K~K~^Tr~^ {K~^ misidenti- 
fied as p) at very low levels consistent with the predictions 
of our detector simulation, but no unexpected structure. 
Prom these studies we calculate that reflections known to 
give broad structures in the vicinity of the yl+ peak, such 
as if~7r'^7r+ (7r+ misidentified as p) contribute 

a number of entries in each bin much smaller than the 
statistical fiuctuations. We also study processes such as 
^ yl+7r+, with the wrong 7r+ included in the A'^ 
candidate, and find their contributions to be negligible. 
In each Xp bin, a linear function describes the mass dis- 
tribution in the data over a wide range away from the 
yl^ peak region, so we use a linear background function 
and perform fits over the range 2235-2335 MeV/c^. 



We first fit the full data sample in each Xp bin in order 

to study mass resolution and bias. These fits yield 
mass values that vary slightly with Xp in a manner con- 
sistent with the simulation and our recent measurement 
of the yl+ mass [30]. The fitted mass resolutions (RMS 
width of the signal function) are shown as a function of 
Xp in Fig. 4. The simulation is consistent with the data 
bin in the data, we fit the weighted invariant mass at low Xp, and is slightly optimistic at high Xp. The effect 



IV. SIGNAL EXTRACTION 



To estimate the number of yl+— >pif 7r+ decays in each 



distribution with a function comprising signal and back- 
ground components. Based on the simulated mass distri- 
butions, we describe the signal with a sum of two Gaus- 
sian functions of common mean value, one of which has 
1.5 times the width and one quarter of the area of the 
other, and correct for a 1.3% residual bias in the fitted 
area. We check the simulated bias by comparing with a 
single Gaussian signal function. The change in the yields 
is 1.2% in both data and simulation. 



of this diff'erence on the efficiency estimate is negligible. 

Next we fix the mean and width of the signal func- 
tion in each Xp bin to values from linear parametriza- 
tions, and perform fits to the on- and off-resonance 

data separately. Dividing the signal yields by the in- 
tegrated luminosity, bin width and branching fraction 
BpK-^+=B{A+^pK-Tr+) = 5.0±1.3% [4] gives the dif- 
ferential production cross sections shown in Fig. 5 with 
statistical errors only. We use the e+e~^cc acceptance 
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FIG. 4: RMS width of the fitted /!+ signal funetion in the 
data (circles) and simulation (squares) as a function of Xp. 
The line represents a linear parametrization of the data. 




FIG. 5: Differential cross sections for A'^ + yl^ production 
in the off- (circles) and on-resonance (open squares) data as 
functions of Xp. The errors are statistical only. 



factor Acc for the off-resonance data and, for purposes 
of this comparison, an average value of A in each bin 
weighted by the relative production rates measured below 
(sec Tables I and V) for the on-resonance data. There are 
two broad peaks in the on-resonance cross section, corre- 
sponding to the contributions from T{4S) decays at low 
Xp and from e+e^^cc events at high Xp. For Xp >0.47, 
the kinematic limit for a B decay including a A'^ and an 
antiproton, the two cross sections are consistent, indicat- 
ing no visible contribution from T(4S') events. 

We extract the cross section for T(4S') decays by re- 
peating the analysis using Ar for both data sets. In 
each Xp bin we then subtract the off-resonance cross sec- 
tion, scaled down by 0.8% to account for the dependence 
of the cross section on the cm. energy, from the on- 
resonance cross section. We divide the off-resonance and 
T{4:S) cross sections by the e+e~^hadrons and effective 
e+e^ T{4:S) cross sections, respectively, to yield the 
differential production rates per event discussed in the 
following sections. 

We consider several sources of systematic uncertainty. 
We propagate the uncertainties on the measured track 
finding and particle identification efiiciencies by recalcu- 
lating s{p, 9) with all efficiencies of a given type varied 
simultaneously and repeating the fits. Tracking gives an 
uncertainty of 2.5% and the particle identification con- 
tributions total 1.5-2.1%, depending on Xp. We obtain 
a 0.9% uncertainty due to the resonant substructure of 
the A'^ ^ pK~'iT'^ decay similarly. The imccrtainty on 
our integrated luminosity is 1.0%. Simulation statistics 
contribute 2-4% where the rate is significantly nonzero. 
Wc check the fitting procedure by floating the signal 
mean and/or width, fixing them to nominal or fitted 
values, using a single Gaussian signal function, using a 
quadratic background function, and varying the bin size. 
All changes in the signal yields are less than the corre- 



sponding statistical errors, and we take the largest change 

in each Xp bin as a systematic uncertainty. Each simu- 
lated bias is varied by ±50%; together they contribute 
0.7% to the systematic uncertainty. The imperfect Xp 
distribution used in the efficiency calculation can afi^ect 
the result if the eSiciency varies over the width of an Xp 
bin. Wc recalculate the efficiency with the input distri- 
bution shifted by plus and minus our bin width to derive 
a conservative limit on any such effect of 0.5-1.9%, de- 
pending on Xp, which we take as a systematic uncertainty. 

We also perform several systematic checks of the re- 
sults. The cross sections measured separately for yl+ and 
which have very different efficiencies at low labora- 
tory momentum, are consistent. Cross sections measured 
in six different regions of cos 6* are consistent with each 
other. Due to the boosted cm. system, these would be 
affected differently by any deficiency in the detector simu- 
lation, especially at low Xp. They also have very different 
Acc values, and these studies indicate that the uncertain- 
ties due to both the production angle model in cc events 
and the value of the boost are negligible. 



V. RESULTS FOR = 10.54 GeV 

A. Baryon Production 

The differential yl+ production rate per hadronic event 
{l/Nqq)[dN'^^/dxp) is tabulated in Table I and compared 
with previous charmed baryon measurements in Fig. 6. 
Wc distinguish between systematic uncertainties that af- 
fect the shape of the cross section and those that affect 
only its normalization. The former include both uncer- 
tainties that are uncorrelated between bins and the parts 
of the correlated uncertainties whose values depend on 
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FIG. 6; Differential ylj production rate per e^e^^qq event 
compared with previous measurements. The error bars in- 
clude statistics and those systematic errors that affect the 
shape. Each experiment has a normalization uncertainty of a 
few percent, and there is an overall 26% uncertainty due to 
the >p/f~7r"'" branching fraction. 



Xp. The uncertainty from the fitting procedure has neg- 
ligible correlation between bins, and those from the par- 
ticle identification and the shift of the simulated distri- 
bution have only very short-range correlations, so wc in- 
clude them in the uncorrelated category. We express the 
error matrix for the efficiency calculation as the sum of 
a diagonal "uncorrelated" matrix and a remainder ma- 
trix, in which the elements of the former are as large as 
possible but no correlation coefficient in the latter ex- 
ceeds unity. The sum in quadrature of the uncorrelated 
uncertainties is listed in the "independent" column of 
Table I. It is typically 3% in the peak region, increasing 
to 10% where the cross section is one-third of its peak 
value, and becoming relatively large at the ends of the 
Xp range. The square roots of the diagonal elements of 
the remainder matrix are listed in the "correlated" col- 
umn of Table I, and included with the independent and 
statistical components in the error bars in the figures. 

All other uncertainties are fully correlated between 
bins and very nearly independent of Xp, so are con- 
sidered experimental normalization uncertainties. They 
total 2.9%, dominated by the track-finding efficiency. 
There is a 26% uncertainty on that also af- 

fects the normalization. The integral of the differential 
rate, taking the correlation in the errors into account, 
gives the total rate, listed at the bottom of Table I 
along with the normalization uncertainties. The prod- 
uct_of the total rate per event and branching fraction of 
Nil ■ ^pK-^+ = 2.84 ± 0.04 (stat.) ± 0.09 (syst.) x 10-^ 
is consistent with, and more precise than, previous mea- 
surements. The normalization uncertainties are not in- 
cluded in any of the figures, and all rates shown assume 
the same value of Bpx--K+- 



TABLE I: At differential production rate per hadronic event 
per unit at 10.54 GcV. The last column includes ex- 
perimental errors that arc correlated between Xp values and 
affect the shape of the distribution. Normalization uncertain- 
ties are given only on the total. 





1 dNf^ 

Nqq dXp 


kj U CLL±0 Li H_vCt± 


Systematic 


Tl CP 

J_ LCtlltll-' 


Error 


Independent Correlated 


0.000-0.025 


0.0007 


0.0009 


0.0004 


0.0001 


0.025-0.050 


0.0033 


0.0016 


0.0012 


0.0002 


0.050-0.075 


0.0008 


0.0030 


0.0011 


0.0001 


0.075-0.100 


-0.0023 


0.0040 


0.0013 


0.0001 


0.100-0.125 


0.0105 


0.0053 


0.0020 


0.0007 


0.125-0.150 


0.0065 


0.0053 


0.0041 


0.0002 


0.150-0.175 


0.0172 


0.0057 


0.0024 


0.0004 


0.175-0.200 


-0.0006 


0.0056 


0.0036 


0.0000 


0.200-0.225 


0.0197 


0.0057 


0.0038 


0.0004 


0.225 0.250 


0.0180 


0.0059 


0.0039 


0.0003 


0.250 0.275 


0.0323 


0.0061 


0.0064 


0.0006 


0.275-0.300 


0.0324 


0.0056 


0.0040 


0.0006 


0.300-0.325 


0.0273 


0.0054 


0.0064 


0.0006 


0.325-0.350 


0.0517 


0.0056 


0.0053 


0.0011 


0.350-0.375 


0.0509 


0.0053 


0.0029 


0.0010 


0.375-0.400 


0.0617 


0.0054 


0.0045 


0.0012 


0.400-0.425 


0.0759 


0.0055 


0.0040 


0.0014 


0.425-0.450 


0.0667 


0.0051 


0.0032 


0.0010 


0.450-0.475 


0.0939 


0.0055 


0.0044 


0.0014 


0.475-0.500 


0.1051 


0.0056 


0.0041 


0.0014 


0.500-0.525 


0.1126 


0.0056 


0.0048 


0.0014 


0.525-0.550 


0.1220 


0.0056 


0.0043 


0.0015 


0.550-0.575 


0.1403 


0.0058 


0.0041 


0.0014 


0.575-0.600 


0.1526 


0.0058 


0.0044 


0.0015 


0.600-0.625 


0.1548 


0.0058 


0.0061 


0.0014 


0.625-0.650 


0.1394 


0.0055 


0.0038 


0.0012 


0.650-0.675 


0.1409 


0.0052 


0.0045 


0.0012 


0.675-0.700 


0.1352 


0.0052 


0.0037 


0.0011 


0.700-0.725 


0.1232 


0.0049 


0.0035 


0.0010 


0.725-0.750 


0.0979 


0.0043 


0.0030 


0.0009 


0.750-0.775 


0.0803 


0.0040 


0.0026 


0.0007 


0.775-0.800 


0.0673 


0.0034 


0.0035 


0.0008 


0.800-0.825 


0.0464 


0.0029 


0.0026 


0.0006 


0.825-0.850 


0.0332 


0.0025 


0.0017 


0.0004 


0.850-0.875 


0.0220 


0.0020 


0.0010 


0.0005 


0.875- 0.900 


0.0161 


0.0017 


0.0017 


0.0006 


0.900 0.925 


0.0079 


0.0013 


0.0007 


0.0003 


0.925 0.950 


0.0049 


0.0012 


0.0006 


0.0008 


0.950-0.975 


0.0018 


0.0006 


0.0003 


0.0004 


0.975-1.000 


0.0003 


0.0001 


0.0002 


0.0001 


Total 


0.0568 


0.0007 


0.0006 


0.0008 


Norm. err. 


0.0016 








BF error 


0.0148 









Assuming the arc produced predominantly in 
e~^e~ — > cc events, the total rate corresponds to a rate 
of Ar^^=0.071±0.003 (exp.)±0.018 (BF) A+ per c-quark 
jet. Roughly 10% of the particles in high-energy jets 
have generally been observed to be baryons [4], and our 
measurement is consistent with 10% of c jets producing 
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FIG. 7: Xp distribution for in e^e^^qq events compared 
with those measured for (a) by BABAR [3] and (b) charmed 
mesons by Belle [2]. 



TABLE II: The first six moments of the /ij Xp distribution 
in hadronic events at v^= 10.54 GeV. 



Stat. Systematic 



Moment 


Value 


Error 


Indep. 


Correl. 


Belle 




0.5738 


0.0061 


0.0049 


0.0032 


0.5824±0.0025 


(4) 


0.3544 


0.0038 


0.0030 


0.0021 


0.3649±0.0034 


(4) 


0.2305 


0.0026 


0.0021 


0.0015 


0.2396±0.0023 


(4) 


0.1560 


0.0020 


0.0015 


0.0011 


0.1630±0.0051 


(4) 


0.1090 


0.0015 


0.0012 


0.0009 


0.1151±0.0020 


(4) 


0.0783 


0.0012 


0.0010 


0.0008 


0.0851±0.0023 



tween different heavy hadrons, and the higher moments 
of the distribution are of theoretical interest. In Table II 
we list values of the first six moments of the Xp distri- 
bution, calculated by summing over bins. They are con- 
sistent with previous measurements from Belle [2]; all 
are 1-2 standard deviations lower, but the moments are 
strongly correlated with each other. The (Xp) value of 
0.574±0.009 is consistent with those measured [2] for 
and mesons, and about 5% lower than those for 
and Z)*+ mesons. 



B. Tests of c-Quark Pragmentation Models 



a c baryon, with a large fraction of these decaying via 
a yl+. All known non-strange charmed baryons decay 
predominantly through a whereas no f}^ states and 
only the heaviest observed Sc states [31] are known to de- 
cay through a so that 70-85% of inclusive charmed 
baryons would be expected to decay through a yl+ in 
current hadronization models. 

The shape of the differential production rate is consis- 
tent with previous results and measured more precisely. 
It is quite hard, as expected, peaking near Xp = 0.6. We 
normalize the rate to unit area to obtain P{xp), and 
compare it with previously measured distributions for 
baryons and D and D* mesons in Fig. 7. The E^ dis- 
tribution is normalized to have the same peak height as 
the yl+ distribution, and, since it was measured on the 
T{4:S) resonance, is shown only above the kinematic limit 
for S-meson decays. The two charmed baryons have sim- 
ilar distributions, with that for the heavier baryon shifted 
up in Xp by roughly 0.05. Although qualitatively similar, 
the D meson distributions show broader peaks than the 
baryon distributions and differ greatly in the way they 
fall toward zero at high Xp. The charmed baryon and 
meson distributions are all much softer than the inclu- 
sive i?-hadron distribution at cm. energies well above bb 
threshold, which peaks around a;p=:0.75 [16]. 

The average Xp value is often used in comparisons be- 



Testing models of heavy-quark fragmentation can be 
problematic since the predictions are usually functions 
of a variable z that is not accessible experimentally, 
such as zi = {E + pii)h/{E + pi\)q, Zi = V\\^Iv\\q or 

^3 = Ph /PHmax{PQ), whcrc p\\ represents a momentum 
projection on the flight direction of the heavy quark be- 
fore it hadronizes. Monte Carlo event generators use sim- 
ilar internal variables, and in some cases can be made 
to produce events according to a given input function 
f{z,(3), where /3 represents the set of model parameters. 
In this way one can test the large-scale features of any 
model, although the detailed structure may not be repro- 
duced exactly. 

We consider the perturbative QCD calculations of 
ColHns & Spiller (CS) [5] and Braaten ct al. (BCFY) [7], 
as well as the phenomenological models of Kartvelishvili 
et al. for mesons (KLP-M) [9] and baryons (KLP-B) [13], 
Bowler [10], Peterson et al. [11], the Lund group [12], the 
UCLA group [14] and the HERWIG group [32]. The lat- 
ter two include heavy quark fragmentation within their 
own generators, and the other seven predict the func- 
tional forms listed in Table III. We implement each of 
these functions within the JETSET generator. 

JETSET uses Zi as its internal variable, but Z2 and z^ 
are very similar at high Xp where we are most sensitive 
to the shape. All distributions are affected by JETSET's 
simulation of hard and soft gluon radiation. 
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TABLE III: Fragmentation models compared with the data. 
Here = ttih + pj_, p± is the component of the hadron 
momentum transverse to the quark momentum, and the Zi 
are defined in the text. 



Model 




Ref. 


CS ([1 


- zi]/zi + e[2 - 


Zl]/[1-Zl]) {1+zf) X 
l/zi - e/\l - zill"^ 


[5] 


BCFY 


22^-1- ^-2; ~ 

[3-3«2-( 3 

+ «i-(i2- 

- zi-{ 9 
+ 34-( 1 


[i — a\Z2 ) X 

- 4d) 

- 23d + 26d2) 
-lld+12d^)(l-d) 

- d+ d^){l-df] 


[7] 


KLP-M 






[9] 


KLP-B 






[13] 


Bowler 




2:3)" cxp(-fem5^/23) 


[10] 


Peterson 


(1A2)(1- 


[l/22]-e/[l-22])-2 


[11] 


Lund 


(iM)(i- 


21)" exp(— 6mi/2i) 


[12] 



We test each model against our measured P{xp) using 
a binned 

n 

_ ^ ^ i^pdata _ pMCyy-l ^pdata _ pMC^ ^^--^ 

where n is the number of bins, Pf"-^'^ {Pf^) is the frac- 
tion of the measured (modelled) distribution in bin i, and 
V is the full error matrix formed from the errors on the 
data (Table I) with the statistical errors on the simula- 
tion added in quadrature to the diagonal elements. For 
each function in Table 111 we minimize this with re- 
spect to the set of parameters /3 by scanning over a wide 
range of possible /? values and generating a large sample 
of qq events at each of several points. We then gener- 
ate additional sets near each minimum. All other model 
parameters are fixed to their default values. 

The functions with one free parameter (CS, BCFY, 
KLP-M, KLP-B and Peterson) all show a single, well 
behaved minimum. The two parameters in the Bowler 
and Lund functions are strongly correlated, and the 
shows a single narrow valley. The UCLA model has in- 
ternal parameters a and b that control production of all 
particles simultaneously; since it has been suggested that 
their values should be different for mesons and baryons, 
we vary them by the same procedure, again finding a 
strong correlation and a narrow x^ valley. HERWIG has 
no free parameter controlling heavy hadron production, 
so we consider only the default parameter values. 

We compare the fitted distributions with the data in 
Fig. 8, and list the x^i the fitted parameter values and 
the average Xp of each fitted distribution in Table IV. 
The parameter values are meaningful only in the context 
of the JETSET (or UCLA) model. The KLP-B, Lund 



TABLE IV: Results of the fragmentation model tests. The 
minimum value, number of degrees of freedom, fitted pa- 
rameter values, and the mean value of the corresponding 
scaled momentum distribution are listed. 



Model 


xVdof 




Parameters 


(Xp) 


JETSET+CS 


227/39 




= 0.135 






0.563 


JETSET-fBCFY 


234/39 


d 


=0.355 






0.560 


JETSET+KLP-M 


219/39 


a 


= 3.05 






0.572 


JETSET+KLP-B 


48/39 


a 


= 7.62 






0.580 


JETSET-hBowler 


52/38 


a 


= 0.93, 


b = 


0.88 


0.583 


JETSET+Peterson 


100/39 


e 


= 0.077 






0.559 


JETSET+Lund 


49/38 


a 


= 1.20, 


b = 


0.71 


0.584 


UCLA 


107/38 


a 


= 2.9, 


6 = 


0.74 


0.584 


HERWIG 


456/40 










0.546 



and Bowler models give the best descriptions of the data, 
with respective confidence levels of 0.15, 0.11 and 0.06. 
However their fitted distributions are systematically be- 
low the data at the lowest Xp values and above the data 
just below the peak region. The UCLA distribution is 
qualitatively similar to these three models but falls more 
rapidly at low Xp, resulting in poor agreement with the 
data. The CS, BCFY and KLP-M models predict dis- 
tributions that are much too broad, and the Peterson 
distribution is also too broad. The HERWIG distribu- 
tion is consistent with the data in the peak region, but 
cuts off too sharply at high Xp. 

The fitted values of the parameter a for the UCLA, 
JETSET+Lund and JETSET-hBowler models arc larger 
than those that describe the production of inclusive light 
hadrons and charmed mesons (a « 1.2 for UCLA and 
0.1 < a < 0.6 for the other two models). Differences be- 
tween baryon and meson distributions have been sug- 
gested on the basis of quark counting [13, 14]. Cross- 
ing of the diagram for leading hadron production in 
e'^e' annihilation gives a deep inelastic scattering dia- 
gram, calculations for which depend on the number of 
spectator quarks, Ng-, in the limit z — >1, one expects 
/(z)(x(l-z)2^'-i [13, 33], and 2iV^-l = 3for baryons. 
This is the form of the KLP-B function, which provides a 
much better description of the data than its counterpart 
for mesons. The UCLA model and the Lund and Bowler 
functions also contain (1 — z)°' terms. For UCLA, the 
fitted value of a = 2.9 is close to 3, as anticipated [14]. 

The models predict P{xp) for primary leading charmed 
hadrons, whereas the data also contain secondary 
charmed hadrons from the splitting of hard gluons, and 
some of the reconstructed are decay products of other 
charmed baryons. Both of these effects are included in 
the JETSET, HERWIG and UCLA models, but it is im- 
portant to consider the effects of possible mismodelling. 
The fraction of qq events containing a gluon splitting 
into a cc pair has been measured at = 92 GeV to be 
about 0.01 [34]. At our lower cm. energy this rate is ex- 
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FIG. 8: The A'^ Xp distribution (dots) compared with the results of the model tests (histograms) described in the text. The 
error bars include simulation statistics. 



pected to be reduced, and the fraction of these that pro- 
duce charmed baryons is expected to be lower than that 
for primary c and c quarks. The JETSET, UCLA and 
HERWIG models predict overall contributions of only 
(0.009±0.004)%, (0.017±0.005)% and (0.034±0.012)%, 
respectively, concentrated at low Xp. The uncertainties 
in our lowest-Xp bins are large enough to accommodate 
such a contribution. Adjusting the models to remove or 
double this contribution does not change the results of 
the fits significantly. 

Currently there are three known Sc states (with 
masses 2455, 2520 and 2880 MeV/c^) that decay into 
yl+TT, and four Ac states (masses 2593, 2625, 2765 and 
2880 MeV/c^) that decay into A+tt+tt". In decays of 
such baryons into a slightly lighter baryon and one or 



two pions, the daughter baryon carries most of the mo- 
mentum, so the effect of any such decay is to soften 
P{xp) slightly without distorting it substantially. In the 
JETSET+Lund, JETSET-hBowler and UCLA models, 
this effect is partially compensated by the fact that the 
heavier baryons are generated with slightly harder dis- 
tributions. Collectively, the effect is to broaden P{xp) 
slightly and shift its average value down. Combining our 
candidates with additional pions in the same event, 
we see clear signals for all of these states, and can com- 
pare the relative contributions to the detected A'^ with 
the simulation. The largest contribution of about 7% 
from i7c(2455) is well simulated in all models, but the 
Z'c(2520) rate is too high by a factor of ^3, and the 
excited Ac states are not in any simulation. Removing 
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two-thirds of the I7c(2520) narrows the simulated distri- 
butions shghtly, but does not improve any value signif- 
icantly. Similarly, adding excited Ac states broadens all 
distributions slightly, with no change in the conclusions 
of the model tests. No Sc or He states are known to de- 
cay to yl+ [4], except two recently reported by Belle [31]. 
The latter are observed at very low rates, so should have 
neghgible effect on P{xp). 



VI. RESULTS FOR r(45) DECAYS 
A. Charmed Baryon Production 

The differential production rate per T(45') decay is shown 
in Fig. 9 and listed in Table V. The errors are as for the 
qq results, with an additional 1.5% normalization uncer- 
tainty due to the qq subtraction procedure. The kine- 
matic limit for T{AS)^ BB—* A^p decays is Xp = 0.47, 
and above this value the rate is consistent with zero. This 
region is omitted from the table and only partly shown in 
the figure. We calculate the total rate by integrating the 
differential rate over the kinematically allowed bins. The 
resulting product of total rate and branching fraction, 



= 4.56 ±0.09 (stat.) ±0.31 (syst.) x 10" 



is consistent with previous measurements [2, 24]. It corre- 
sponds to 7VX; = 0.091 ±0.006 (exp.)±0.024 (BF) A+ per 

T{AS) decay, and 4.5±1.2% of S^/B" decays including 
a yl+ baryon, assuming the T(4S') decays predominantly 
to BB. 

Our results on the shape are consistent with, and more 
precise than, previous results, which are also shown in 
Fig. 9. The Xp distribution is quite soft. In particular, the 
data drop rapidly above the peak and are consistent with 
zero above Xp^ 0.35. This is the range expected for quasi- 
two-body decays into a yl+ or Ec plus an antibaryon such 
as a n or A, and includes much of the range expected 
for decays involving one or two additional pions. The 
B^^ A^p decay has been observed [21] at a very low rate 
consistent with our inclusive data. 

A soft Xp distribution was also seen in our recent study 
of baryons [3], but the statistics of that study did 
not allow a meaningful direct subtraction of the con- 
tribution from qq events. As an exercise, we assume a 
smooth distribution from qq events by choosing an em- 
pirical function that describes both the A'^ data in Sec. V 
and the high-Xp data (Fig. 7a). We fit this function 
to the high-Xp data and subtract the result in all Xp 
bins. The resulting approximate Xp distribution for 
in T(AS) decays is also shown in Fig. 9, normalized to 
have roughly the same peak height as the vl+ data. It 
is also quite soft, similar in shape to the yl+ but shifted 
slightly downward in Xp, and consistent with zero above 
Xp« 0.35. Because of the e+e"^ cc subtraction proce- 
dure, the error bars cannot be compared with those for 
the A'^ data, but the noted features do not depend on 
the details of the subtraction procedure. 



• A^ BaBar 

• A^ Belle 
o A^ CLEO 




FIG. 9: Differential yl+ production rate per T(4S') decay com- 
pared with previous measurements. Normalization errors are 
not shown, and the S'^ rate is normalized to match the peak 
At rate. 



B. Model Tests 



Existing models of B meson decays into c baryons were 
developed with little data on the Xp distribution. In 
Fig. 10a we compare P{xp) (normalized over the range 
Xp < 0.475) with the predictions of the models in the 
JETSET generator [28] and our internal generator [29]. 
The former has been tuned to measured multiplicities 
and momentum distributions of stable B decay products, 
and predicts a distribution that is much too soft. The 
latter includes many measured exclusive decays involv- 
ing D mesons and postulates analogous few-body decays 
involving charmed baryons. Its predicted Xp distribution 
is similar in shape to the data, but shifted to higher Xp 
values. Neither generator produces simple (quasi-)two- 
body decays such as B-^A^p, A^n, ^ EcP- 

In Fig. 10b we compare the same data with simulated 
events of the type B^A^p{nnr) for selected values of m, 
the number of pions in the decay in addition to the Af 
and antiproton. The distributions arc insensitive to the 
charges of the pions or B meson, or to replacing the an- 
tiproton with an antineutron. Decays via a A or strange 
antibaryon are not included; they give A^ distributions 
only slightly different from those shown with m — 1 pi- 
ons. For m= 1, 2, 4 and 6, the distributions shown are 
from the JETSET simulation; phase space decays give 
similar distributions. The spread in the distribution for 
m = is due to the finite momentum of the B meson in 
the e+e" cm. frame. The measured P{Xp) is described 
adequately by the simulation with m = 4. Adding con- 
tributions from 777, = 3 and m = 5 improves the of ^ 
comparison with the data, but no further contributions 
are helpful. That m is restricted to such a narrow range 
suggests that different types of decay modes are needed. 

An intriguing possibility is that there is a large con- 
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TABLE V: At differential production rate per T{4S) decay 
per unit scaled momentum (to \/s = 10.58 GeV). The last 
column includes those experimental errors that are correlated 
between Xp values and affect the shape of the distribution. 
Normalization uncertainties are given only on the total. 



Xp 

Range 


r Ac 


Statistical 


Systematic 


Error 


Independent Correlated 




0.000-0.025 


0.0123 


0.0023 


0.0040 


0.0013 


0.025-0.050 


0.0563 


0.0073 


0.0135 


0.0046 


0.050-0.075 


0.1963 


0.0110 


0.0443 


0.0182 


0.075-0.100 


0.3317 


0.0152 


0.0194 


0.0215 


0.100-0.125 


0.3686 


0.0193 


0.0196 


0.0294 


0.125-0.150 


0.4555 


0.0208 


0.0285 


0.0205 


0.150-0.175 


0.4669 


0.0221 


0.0238 


0.0167 


0.175-0.200 


0.4482 


0.0201 


0.0240 


0.0151 


0.200-0.225 


0.3863 


0.0211 


0.0454 


0.0132 


0.225-0.250 


0.3372 


0.0204 


0.0339 


0.0109 


0.250-0.275 


0.2315 


0.0197 


0.0246 


0.0084 


0.275-0.300 


0.1447 


0.0183 


0.0327 


0.0064 


0.300-0.325 


0.1255 


0.0171 


0.0227 


0.0055 


0.325-0.350 


0.0365 


0.0175 


0.0171 


0.0041 


0.350-0.375 


0.0336 


0.0160 


0.0072 


0.0037 


0.375-0.400 


0.0031 


0.0163 


0.0155 


0.0043 


0.400-0.425 


-0.0117 


0.0163 


0.0080 


0.0042 


0.425-0.450 


0.0268 


0.0150 


0.0081 


0.0044 


0.450-0.475 


-0.0096 


0.0158 


0.0097 


0.0048 


Total 


0.0910 


0.0019 


0.0026 


0.0049 


Norm, error 


0.0029 








BF error 


0.0237 









tribution from decays involving both a c baryon and 
an anti-c baryon. The decays 5° ^ — [22] and 
B~ — > A'tA~K~ [23] have recently been observed, and 
we have previously measured an unexpectedly high rate 
of inclusive A~ in B~ decays [20]. As an exercise, 
wc model two-body two-c-baryon decays based on the 
simplest internal W diagrams, i.e. those of the forms 

and 5° 5^{I!^)IJ^, where the parentheses indicate 
Cabibbo-suppressed modes. Here, Sc represents any of 
the states ='^(2470), S'^(2570) or S'c(2645), Uc represents 
i;c(2455) or i;c(2520), Ac represents ylc(2285), ylc(2593) 
or ylc(2625), and we consider all kinematically allowed 
combinations at relative rates determined by phase space. 
We decay all Ec baryons into A'^n and all excited Ac 
baryons into vI+tttt, with 73% of the vl+(2593) baryons 
decaying through a ScTt intermediate state and all others 
via phase space. 

The Xp distribution of the yl+ from this simulation 
is compared with the data in Fig. 10c. Although it is 
too narrow to describe the data completely, it appears 
that such processes could contribute substantially to the 
overall rate. Combining this simulation with those for 
the yl+|5(m7r) modes and assuming a smooth, broad 
distribution of m, we can describe the data with as much 




FIG. 10: P{xp) for At from T(4S) decays (dots) compared 
with the simulations (lines) described in the text. The data 
are normalized to unit area, as are the BABAR simulation in 
(a) , the B — > Atp4n simulation in (b) and the composite sim- 
ulation in (c); all other simulations are normalized arbitrarily. 



as a 50% contribution from these two-c-baryon decays. 
For example, the "composite model" in Fig. 10c com- 
prises 35% two-c-baryon decays and (12, 25, 12, 9, 7)% 
of m=(2, 3, 4, 5, 6), and provides an excellent descrip- 
tion of the data. Decays with two charmed baryons and 
additional pions or kaons could also contribute at low Xp, 
shifting the m distribution downward. Measurements of 
many exclusive baryonic B decays, including both one- 
and two-c-baryon modes, are needed to understand the 
dynamics in detail. 
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VII. SUMMARY AND CONCLUSIONS 

We use the excellent tracking and particle identifica- 
tion capabilities of the BaBAR detector to reconstruct 
large, clean samples of /1+ baryons over the full kine- 
matic range. We measure their total production rates 
and inclusive scaled momentum distributions in both 
e+e" hadrons events at ^/s = 10.54 GeV and 1^(45) 
decays. Our results are consistent with those published 
previously and more precise. 

In e+e~^gg events we measure a total rate per event 
times branching fraction into the pK~n~^ mode of 

Nfc ■ ^pK-n+ = 2.84 ± 0.04 (stat.) ± 0.09 (syst.) x 10"^, 

where the first error is statistical and the second sys- 
tematic. The imccrtainty on the total rate per qq event, 
Ar9|= 0.057 ± 0.002 (exp.) ± 0.015 (BF), is dominated by 
the uncertainty on Bpx-Tr+ ■ The corresponding value of 
= 0.071±0.003 (cxp.)±0.018 (BF) A+ per c-quark 
jet is consistent with the hypothesis that roughly 10% of 
c jets produce a c baryon and a large fraction of these 
decay via a A'^. 

The scaled momentum distribution peaks at XpWO.6. 
It is similar in shape to those measured previously for 
D and D* mesons, but peaks more sharply and drops 
toward zero more rapidly as 1. We measure an 

average value for yl+ of 

{xp) = 0.574 ± 0.006 (stat.) ± 0.006 (syst.), 

which is consistent with values measured for ground state 
D mesons, but about 5% lower than those for D* mesons. 
We use this distribution to test several models of heavy 
quark fragmentation, none of which provides a complete 
description of the data. The baryon-spccific model of 
Kartvelishvili et al. and the models of Lund and Bowler 
have acceptable values, but all show a steeper slope on 
the low side of the peak than the data. The UCLA model 
shows similar qualitative features, but worse agreement 
with the data. The HERWIG model is far too narrow, 
and all others are too broad. In previous model tests 
using specific c mesons [2] and inclusive b hadrons [16] 
(a mix of roughly 90% mesons and 10% baryons), the 
Lund, Bowler and Kartvelishvili models generally gave 
the best description of the data, and UCLA described 
the 6-hadron data, whereas the other models showed dis- 
crepancies similar in form to those reported here. The 
Kartvelishvili and UCLA models postulate different spec- 
tra for mesons and baryons. Their strong preference for 



their respective baryonic forms, combined with the ob- 
served differences in shape between our yl+ spectrum and 
previously measured D meson spectra, indicate a differ- 
ence in the underlying dynamics. 

In T{AS) decays, we measure a total rate per event 
times branching fraction of 

■ l3pK-n+ = 4.56 ± 0.09 (stat.) ± 0.31 (syst.) x 10"^, 

corresponding to N^^ = 0M5 ± 0.003 (exp.)±0.012 (BF) 
per B°/B~ decay. The spectrum is softer than predicted 
by our B decay model, and much harder than that pre- 
dicted by JETSET. It can be described by models of 
B— > A^p{rmr) decays only if the m = 3-5 contributions 
dominate. Alternatively, a model including a large con- 
tribution from decays involving both a charmed and an 
anti-charmed baryon can describe the data in conjunc- 
tion with a broad distribution of m. Additional studies 
of exclusive modes are needed to understand the details 
of B-meson decays into baryons. 
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